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Abstract 

We study the quantum foundations of a theory of large amplitude collective motion for a 
Hamiltonian expressed in terms of canonical variables. In previous work the separation into 
slow and fast (collective and non-collective) variables was carried out without the explicit 
intervention of the Born Oppenheimer approach. The addition of the Born Oppenheimer 
assumption not only provides support for the results found previously in leading approxima- 
tion, but also facilitates an extension of the theory to include an approximate description 
of the fast variables and their interaction with the slow ones. Among other corrections, one 
encounters the Berry vector and scalar potential. The formalism is illustrated with the aid 
of some simple examples, where the potentials in question are actually evaluated and where 
the accuracy of the Born Oppenheimer approximation is tested. Variational formulations of 
both Hamiltonian and Lagrangian type are described for the equations of motion for the slow 
variables. 



1 Introduction 

For the past decade, the authors and their collaborators have been involved in the 
development of a theory of large amplitude collective motion. This enterprise began 
with a pair of papers |], on possible quantum foundations for such a program. Soon 
thereafter || it was realized that the leading approximation, to which most previous 
studies had been confined, was classical in nature, (with subsequent requantization) 
and, as a consequence, almost all later systematic theoretical development, reviewed 
in Ref. ||, was based on the study of this limit. 

Quantum corrections can be important in selected circumstances, however. Con- 
sequently, the elements of a systematic method for including quantum corrections by 
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expansion about the classical limit has recently been developed 0, distinct from the 
approaches found in our initial papers. We consider it worthwhile, nevertheless, in the 
light of subsequent developments, especially within the context of the Berry phase idea 
[@, 0, to update, augment and improve the work described in those initial efforts and 
to tie them, where possible, to more recent work. In a paper currently in preparation, 
we shall contrast the nuclear-physics foundations of the method described in this paper 
with that found in Ref. |5[], both approaches requiring elaboration compared to what 
has already been published. 

In Sec. H we give a condensed but at the same time more precise account of the 
essential theoretical content of Ref. namely, a method for separating and identi- 
fying a collective subspace for a specified class of Hamiltonians. This is done in an 
approximation that suppresses any coupling to non-collective degrees of freedom. As 
opposed to the corresponding problem in molecular physics, the Hamiltonian is gener- 
ally given in a form where the coordinates are not appropriate to make the separation 
of the Hamiltonian into a collective and a non-collective part. The discovery of new 
coordinates that will effect the separation is an essential part of the problem. In the 
regime of large amplitude collective motion, a reasonably complete theory that can be 
implemented f|, || has been developed so far only for the case that the transformations 
to new coordinates and momenta are restricted to point transformations. 

In Sec. [|, we describe how the previous theory may be extended to include the 
interaction between the fast and slow variables. This is done by means of a standard 
Born-Oppenheimer (BO) representation of the states of the slow variables. Though 
some of the ensuing details are similar to those encountered in the molecular problem, 
leading, for example, to the occurrence of Berry potentials || [7], others are charac- 
teristic of the problem of large amplitude collective motion. In order to evaluate the 
corrections found to the potential energy of the collective variables, it is necessary to 
obtain wave functions for the fast variables. This is done by an extension of the anal- 
ysis carried out in Sec. 0. At the same time this development provides a justification 
for some of the assumptions of the previous analysis. 

In Sec. |], we present material that has no direct counterpart in any previous work 
in this field, outside our own. One of the results of Sec. § was to establish a complete 
effective quantum mechanics in the collective subspace, i. e., we not only computed 
an effective Hamiltonian operator in the collective subspace, but also proved that it 
was expressed in terms of canonical variables. In Sec. |] we study variational principles 
for the associated Heisenberg matrix mechanics. If we put aside, momentarily, special 
difficulties associated with the curved spaces that we allow in our general formulation, 
we find quantum analogues of variational principles both for Hamilton's and for La- 
grange's equations. The connection between the two forms is the standard one, even 
when Berry phase terms are included, in contrast to a result based on the use of path 
integrals |7]]. For the general case of a curved space the situation is more involved, and 
only a variational principle for Hamilton's equations is presented and analyzed. 
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The corrections to the potential energy of collective motion found in Sec. |3] are 
of three types. There is first a "quantum" potential arising from the curved space 
(position-dependent mass tensor). This term cannot be evaluated fully with the tools 
developed in our previous work || . In practical applications to systems with N degrees 
of freedom, this term is of order iV~ 2 , and therefore should be small for a many-body 
system. A second term consists of the eigenvalue for the state of the fast variables for 
a fixed value of the slow variables. In our approximation this state is that of a set 
of independent harmonic oscillators, with frequencies dependent on the instantaneous 
value of the collective coordinates. We are mainly (though not exclusively) interested 
in the zero-point motion associated with the ground state. This contribution, which is 
only of order N^ 1 compared to the leading term has been studied previously in simple 
examples ||. A third correction, referred to as the Berry scalar potential, occurs 
together with the Berry vector potential, both appearing in the collective Hamiltonian 
as a consequence of the BO approximation. 

Examples in the literature relevant to the approach described in this paper are not 
numerous. The examples studied by Bulgac || |10| have been stimulating, but a bit 
too simple for our purposes. On the other hand, the work by Girard, LeTourneux, 
and Vinet jjllj] on the dipole-quadrupole problem in spherical nuclei is too complex 
to satisfy our initial needs. Sec. || is devoted to the study of several "in-between" 
models sufficiently simple that the BO approximation can be carried out, including 
the calculation of the Berry potentials, and compared with more exact results. Some 
technical details associated with the material of Sees. || and [3], respectively, are provided 
in two appendices. 



2 Summary of previous results 

In this paper we base our study on a Hamiltonian, that with the help of the summation 
convention, takes the form 

H = H(C, 7r a ) = ^K, + 7(0, (2.1) 

that describes N coordinate and momentum pairs, £ a and 7r a , a = 1...N, that satisfy 
canonical commutation relations, (h — 1), 

[?,<*!,]= i6%. (2.2) 

We thus allow for a curved space described by a (reciprocal) mass tensor B a/3 . Not 
only is this formulation of interest for a range of applications outside of nuclear physics 
[§, H, but it is also possible |12| to transform the large amplitude problem of nuclear 
collective motion into this form as the basis for further development. Associated with 
this Hamiltonian is a form for the scalar product in Hilbert space that is discussed in 
appendix A. 
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The problem of immediate interest is to find a decomposition of the operator H 
into two parts, exactly in rare cases, approximately under most circumstances, that 
describe slow and fast (or collective and non-collective) degrees of freedom. A basic 
difficulty of the class of problems that interest us is that we cannot assume that this 
separation occurs for the initial choice of coordinates. Instead, we assume that it is 
possible to effect such a decomposition, at least approximately, by means of a point 
transformation of the form 

r = g a (x), 

x = {Q\q a }, 
i = 1...K, 

a = K + 1...N, (2.3) 

that is locally invertible, 

= (2-4) 

Though we have discussed, in the past |3|, [|], the possible interest of going beyond 
point transformations in the study of large amplitude collective motion, it is only for 
this simpler class of transformations that a substantial theoretical underpinning can 
be claimed to exist, a foundation that we wish to widen in the present work. (Of 
course, there is a substantial literature based on transformations that are polynomial 
in coordinates and momenta that are applicable to anharmonic vibrations (see, for 



instance [13], [bp), but this is not the subject of the current paper.) 

A first step is to carry out a formal transformation of the Hamiltonian, H, to the 
new variables. For the potential energy, we have, trivially enough, 

V(0 = V(g a (Q,q)) = V(Q,q). (2.5) 

Turning then to the consideration of the kinetic energy, T, 

T= ^,{7^,^(0}}, (2-6) 

this expression can be transformed to the new variables with the aid of the relations 
I 

*■« = \{f%P& (2-7) 
B v> = (2.8) 



The derivation of ( |2.7| ) is given in appendix A, whereas ( |2.8| ) confirms the tensorial 
character of the mass tensor. In terms of the new variables, the kinetic energy conse- 
quently takes the form 

T = l{p„ {Pu, B^(Q, q)}} + U(Q, q), (2.9) 
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where the second term, which is specifically a quantum potential arising from the 
non-commutativity of coordinates and momenta, has the form 



8U(Q, q) = [f^9l9%B^X - [f% l9 lU%B vX . (2.10) 

For further work, the momenta are also divided into collective and non- collective sub- 
sets, 

P» = {Pi,Pa}. (2-11) 

We now remind ourselves that the initial, even primary, goal of the present consider- 
ations is to identify a piece of the transformed Hamiltonian as collective. This collective 
Hamiltonian, H c , should depend only on the variables Q\ Pj, that we suppose to be 
canonical pairs. The simplest possible choice would appear to be the restriction of the 
full Hamiltonian to the values q a = 0, p a = 0, and thus (see below) 

H c = ^{P u {Pj, B*(Q, 0)}} + V(Q, 0). (2.12) 

This choice is purely formal until we specify a procedure for determining the unknown 
functions V(Q,0) = V(Q) and B lj (Q,0) = B lj (Q) in terms of the elements of the 
original Hamiltonian and of the point transformation ( |2.3| ) and ( |2.4| ). The criteria for 
such a determination must be chosen so as to somehow minimize the coupling between 
the collective and non-collective variables. 

The simplest approach to this problem is to consider Hq, ( |2.12|) , as the first term 



of a Taylor expansion in powers of the non- collective variables, g a , remembering at the 
same time that we have a polynomial of degree two in the non-collective momentum 
operators, p a . Assuming that higher order terms are successively of lesser importance, 
basically the same assumption that governs the treatment of all variables in the small 
vibrations domain, we adopt the same assumption that defines optimum decoupling 
in the classical limit 0], namely that the terms linear in the non-collective variables 
vanish. This yields the two sets of conditions 

V a {Q) + l -{P t ,{P 3 ,B%Q)}} = 0, (2.13) 

B at (Q) = 0, (2.14) 

recognizable as the coefficients of q a and p a , respectively. The operator character of 
these conditions is apparent only from the second term of ( ^.13[ ). However, if we intro- 
duce the Wigner transform of these equations, because of the special symmetrization 
that we have adopted for the kinetic energy, they reduce exactly to the classical condi- 
tions studied extensively in our previous work There, starting from these classical 
conditions, we have described and illustrated suitable algorithms for the determination 
of the elements of the collective Hamiltonian. 
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With the determination of the elements of the Hamiltonian, ( 2.12|) , the theory is 
completed by the assumption that the eigenstates of this operator are defined on a 
K-dimensional Hilbert space of (round-bracketed) states |n), 



n)= J dQ\Q){Q\n), 



(2.15) 



and that these states accurately model a corresponding subspace, \n), of states of the 
full Hilbert space. The theory described above is very closely the theory that has been 
applied in a number of early applications, such as [[15], [TB|. A major aim of the present 
work is to critique and generalize this procedure. 

We emphasize that in this program, the elements of the collective Hamiltonian are 
determined from a K-dimensional manifold in configuration space, parametrized by 
the coordinates {Q}, whereas the quantum mechanics is determined in a Hilbert space 
defined by the same coordinates. Upon further reflection, it may strike the reader 
that this is a strange result, or, at least, a very limited one. The reason is that the 
true collective states, \n), though they define a subspace of the full Hilbert space, 
are nevertheless states in this space and therefore depend also on the non-collective 
coordinates, at the very least through the zero-point motion of the latter. Though 
we have been able to neglect this dependence in the approximation considered so far, 
this will no longer be true if we seek to include the effect of the coupling to the fast 
variables on the properties of the collective states, or in directly describing the motion 
of the fast variables. 

3 Generalization of the formalism to include the 
effects of the fast variables 

In this section, we set ourselves two tasks. We shall first find the forms of the leading 
corrections to the collective Hamiltonian given above, arising from the coupling of the 
fast to the slow variables. This will be done with the help of the standard Born- 
Oppenheimer (BO) approximation. In order to evaluate these corrections explicitly, 
however, we shall then commit ourselves to further approximations for the dependence 
of the full Hamiltonian on the fast variables, that yield a normal mode description of 
the latter. 

To generalize the restricted ideas of the previous section, we thus introduce the 
BO picture into the theory of large amplitude collective motion. As stated above 
this amounts to an extension of the complexity of structure allowed for the states of 
the collective subspace so as to take into account the influence of the fast variables. 
We first assume, more generally, that the collective states |n) have coordinate space 
representatives 



(QY\n) 



{Q,q\n) 
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= £(QM[gk:Q], 

V 



(3.1) 



where the index v is not to be confused with its previous use as a coordinate index. 
Here and for the remainder of this paper, we adopt a notation where angular brackets 
indicate states in the full Hilbert space, square brackets states in the space of fast 
variables (though dependent parametrically on the slow variables, as denoted by : Q 
in the state vector), and parentheses states in the collective space. We suppose that 
for fixed Q, the states [g|z/:Q] are a complete set of functions for the fast coordinates, 

Y,[q\v.Q][v:Q\q>]=5{q-q'). (3.2) 

Though for the moment we have not specified the equation of which they are the 
solutions, we shall be able to do so, at least approximately, as a consequence of the 
developments to be carried out in this section. 

For the remainder of the current discussion, we shall consider the simplest case 
in which the fast variables occupy, for any given value of the slow variables, their 
state of lowest energy, which is assumed to be non-degenerate. For the case to be 
studied here, v takes a single value denoted by zero. The coefficient function, usually 
denoted in this case by (Q\n), that appears in ( p.l|) , can then be identified as the 
wave function for the collective motion, and this identification will agree with the one 
that has been made in the previous section, where the description of large amplitude 
collective motion was not tied to the BO approximation. There we emphasized a 
connection between the decoupled Hilbert space of the collective coordinates and the 
crucial existence of a i^-dimensional decoupled coordinate manifold described by the 
functions g a (Q,q = 0). The approximate separability of the associated Hilbert space 
is based on the assumption that the functions describing the fast variables are 

confined to a narrow region in q space in the neighborhood of the collective surface. If 
this picture holds, we may expect the mathematical details to work out reasonably. 

Adopting the BO approximation, we set ourselves the task of finding an effective 
Hamiltonian to describe the motion of the collective variables. This operator is defined, 
though not yet operationally, by means of the equation, 

(n'\H eS (Q,P)\n) = (n'\H^,n)\n). (3.3) 

This relation will assume the status of a definition of the collective Hamiltonian within 
the space of slow variables and if e ff recognized as a generalization of He only after we 
specify how to eliminate the fast variables from the right hand side. The procedure 
that we shall follow is closely akin to the traditional BO approach, with characteristic 
differences arising from the facts that at the beginning we cannot specify which are the 
slow and which the fast variables, and that the treatment of the fast variables comes 
as a kind of afterthought, dependent in detail on the prior treatment of the collective 
variables. 
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Let us start with the potential energy, 



V(0 = V(g a (Q,q)) = V(Q,q), (3.4) 
and evaluate the associated piece of (|3.3|) 

(n'\V(0\n) = (ri\V eff \n) = [ dQdq(n'\Q)[0:Q\q]V(g a (Q,q))[q\0:Q](Q\n), (3.5) 



within the BO approximation. The only feasible way, in general, of integrating out the 
fast variables is to expand V in powers of q a , 



V{Q, q) = V{Q) + V a q a + \V ab q a q b + ... , (3.6) 

leading to 

V eff (Q) = V(Q) + V«(Q) + V^(Q) + ... , (3.7) 
where assuming that the function [g|0:(5] is normalized, 

V(Q) = V(Q,0), (3.8) 
V (1 \Q) = V a J\[q\0:Q]U 

= V a (Q)(q a ) Q , (3.9) 
V^(Q) = lV ab (Q)(q a q b ) Q . (3.10) 

We thus see that the leading term is independent of the wave function for the 
fast variables, coinciding with the standard result for the potential energy of large 
amplitude collective motion |], |J presented in the previous section. The computation 
only requires the form of the collective submanifold, £ Q = g a (Q,0), which can be 
determined by well-defined procedures (§]. To go beyond this lowest order, we need, 
besides the wave function of the fast variables, to rearrange the expansion of V, as 
explained in appendix B. There it is shown that if we wish to interpret the small 
quantities q a as components of a vector, we must replace the ordinary second derivative 
by a covariant second derivative, V.^, that can be computed from known or calculable 
quantities according to the equations 

V a (Q) = V a g% (3.11) 
V ab (Q) = V aP g a a g p b (3.12) 

These equations are simplified by introducing one of the decoupling conditions 
that follow from Eq. ( p. 13 ), namely that ( 3 . 1 1| ) should be zero. (Here we are assuming 
that the two terms of ( |2.13| ) vanish separately. Conditions for this to be true and 



modifications necessary when it is not have been discussed in Ref. 0]) The third term 
V^(Q), is the first, A^V(Q), of a sequence of contributions that we shall identify as 
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the leading corrections to the potential energy. Further discussion of the evaluation 
of this term and of the additional correction terms, to be identified below, will be 
continued later in this section, following the identification of all the pieces. 

The remaining terms will arise from the study of the kinetic energy, given after 
transformation of coordinates by Eqs. ( |2.9j ) and ( |2.10| ). In order to integrate out the 
fast variables in the contribution that these terms make to H e ff, consider the first term 
of the kinetic energy. We expand the mass tensor in powers of q, and keep initially 
only the leading term B^ U (Q, 0) = B^ U (Q). In this approximation, we first restrict the 
study to the contribution of those terms where the indices /i, v take on values i, j in 
the collective subset. 

As a preliminary to this calculation, we study the simpler object 

(n'KPOeffln) = (n'\Pi\n) 

= (n'lAW- (3.13) 
Here Pj is the collective momentum operator identified in the previous section, and 

Ai = iJdq[0:Q\q]di[q\0:Q], (3.14) 

where <9j means partial derivative with respect to Q l . Notice, however, that if we 
calculate straightforwardly, 

(n'KPiPjUln) = (ri\PiPj\n) 

= (n'KPiPj-AiPj-AjPi + SiJln), (3.15) 

where 

Sij = -[ dq[0:Q\q}(d 2 /dQ l dQ^[q\0:Q}. (3.16) 

As a consequence of this result, we are reassured that the effective value of zero is 
zero, i. e., 

= (n'\[P l ,P J }\n) = (n'\[P l ,P J }\n). (3.17) 

It is easy to see that the remaining canonical commutators also project without change. 

It is useful to rewrite (|3.15|) in a form that makes contact with standard results. 
We have 

(ic^M^Qli = Y,W-Q\^)iQ))^ (3-i8) 

(Aj) u > u = i J <k[i/:Q\q\dj[q\w.Q\. (3.19) 
Because of the Q dependence of the matrix element, it now follows that 

= (id^Aj + AiAj + S' ip (3.20) 
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where 



(^)oo = A 3 , (3.21) 
5« = E(^W ( 3 - 22 ) 

Consequently, we may also rewrite Eq. (|3.15|) as 

(n'KP^eflk) = (n'\(D i D j + W- (3-23) 

In the simplest case, where [0 : Q\q] is a real wave function, and it follows that Aj 
vanishes, the contribution S'^ remains to be taken into account. 

We are now in a position to apply to the computation of the collective kinetic 
energy the same reasoning as just carried out for a product of momentum operators. 
Making use of the analogue of ( |3.23| ), the result is 

±{D i ,{D j ,BU(Q)}} + kS' ij B«. (3.24) 

where the second term can be incorporated into the collective potential energy as as a 
second such contribution, A^V(Q). 

An additional contribution of this type is obtained by setting q a = in Eq. (|2.10 ), 

A^V(Q) = U(Q,0), (3.25) 

where U(Q,q) is the quantum potential defined in Eqs. ( |2.9| ) and ( |2.1(J| ). 

It remains for us to discuss the contributions from T that depend on B m (Q) that 
"mix" the collective and non-collective indices and those that depend on the non- 
collective mass tensor B ab (Q). The former can be neglected because of one of the 
decoupling conditions, Eq. ( |2.14j ). The leading contribution of the latter is seen to be 
another contribution to the potential energy, 

A^V = \B a \Q)(p aPb ) Q . (3.26) 

To summarize our findings, we have derived the following effective Hamiltonian, 

H ea = {D j: B ij (Q)}} + V(Q) + AV(Q), (3.27) 

where AV(Q) is the sum of four terms that summarize the leading quantum corrections 
including the coupling to the fast variables, 

4 

AV = Y / ^ (i) V, (3.28) 

i=l 
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given respectively in or in relation to Eqs. ( |3.10| ), ( p.24[ ), ( |3.25| ), and ( |3.26| ). Let 
us contrast this result with the corresponding form appropriate to the more familiar 
molecular case. The introduction of a curved metric aside, the main difference is that 
in the molecular case, the ground state wave function of the fast variables, for a fixed 
value of Q, may be assumed known, and its eigenvalue, eo(Q), together with A^V, 
contained in Eq. ( |3.24|) constitutes the collective potential energy || [7j . 



In the present instance we cannot assume that we know the Hamiltonian of the 
fast variables, except in an approximate sense that we now discuss. We need to extend 
the considerations of the previous section that led to the definition of the quantum 
collective Hamiltonian operator He- Instead of specializing the transformed Hamilto- 
nian operator to the values q a = p a = 0, we now retain terms up to second order in 
these variables. In this treatment, we may replace U(Q, q) by U(Q, 0), since this term 
is already a small correction, and the remaining terms linear in the fast variables may 
be dropped because of the decoupling conditions. To the specified accuracy, we obtain 
the following quantum Hamiltonian, 

H = H C + U(Q) + H NC , (3.29) 
H NC = \paP h B a \Q) + \q a q b V ab (Q). (3.30) 

In appendix B, we point out that once the elements of He have been determined, 
one can choose a set of fast variables and calculate the matrices that appear in (|3.30|) . 
For each value of Q, the Hamiltonian H^c represents a standard normal mode problem. 
Let c a , cj, be normal mode destruction and creation operators, Q a the corresponding 
frequencies, and h a = c^ a c a . Assuming local stability, i. e., fl a real and positive, we 
have 

Hnc{Q) = + \)^{Q). (3.31) 

a 

The practical importance of the contribution of (|3.31|) has been noted in several appli- 
cations carried out in the past g, [T5|j . 

The quantum Hamiltonian, expressed in terms of the optimum choice of variables, 
but in a restricted approximation, has thus been found. Except for the term U(Q), it 
has been expressed in terms of elements that can be calculated. Evaluation of U(Q) 
appears to require some properties of the point transformation that have not so far been 
studied. Comparison of approximate calculations (with this term omitted) with exact 
calculations indicates that it is probably a small correction. As previously remarked, 
it is of order N~ 2 compared to the main term in the potential energy, where N is the 
number of degrees of freedom. We shall make no further allusion to this term in the 
present paper. 

Let us now return to a discussion of the correction terms, AV by which H e s, Eq. 
( |3~27|) , differs from H c . We have just discarded A (3) K The sum A {1) V + A^V has 



been seen in (|3.31| ) to be a sum of oscillator terms in the approximation considered, 



its contribution then depending, naturally, on the state of motion of the fast variables. 
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As stated above, in this section we shall consider only the lowest energy state for these 
variables, so that the contribution of this term is just the zero-point energy. 

It remains for us to discuss the contribution of the term A^ 2 ^, the Berry scalar 
potential, associated with the projection of the kinetic energy onto the collective sub- 
space. This contribution goes together with that from the "vector potential", A4. In 
Sec. ||, we shall study examples where these terms contribute. In fact, the manner 
in which both the vector and scalar potentials contribute is best studied within the 
context of these illustrations. 

In this section we have thus gone as far as we shall in defining a quantum theory 
of the slow variables, having derived an effective Hamiltonian and shown that it is 
expressed in terms of canonical variables. Thus we are also free to study this problem 
with the help of the Heisenberg equations of motion. In the next section we shall 
consider variational formulations of these equations. 

4 Variational principles; Lagrangian formulation 

Having defined an effective quantum theory within the collective subspace, we can 
imagine that we study this theory by means of Heisenberg's equations of motion, 

-z[Q\#eff] = J^efr, (4.1) 
i[P t ,H«) = (4.2) 

The purpose of this section is to show that these equations can be derived from a 
variational principle, the so-called trace variational principle, involving H e g, or an 
operator closely related to it; it is to show, further, that a reworking of this principle 
leads to a second version of the variational principle, involving an effective Lagrangian, 
L e g, that is defined in the natural way as the Legendre transform of H s r. According to 
Moody, Schapere, and Wilzek JF], when L e s is defined by means of a path integral, it is 
not exactly the Legendre transform of H s r. We shall not encounter any such difficulty 
as long as the mass tensor is independent of coordinates, the only case studied by 
previous authors. When the mass tensor enters in full generality, the formulation of a 
variational principle of the type we have in mind is not quite as straightforward as in 



previous instances that we have studied [17 



In order not to have to deal with all subtleties at once, let us first consider the case 
that the mass tensor does not depend on coordinates. In that case we may utilize the 
standard version 

5Tr{H cS - iA[Q\ P,}} = 0. (4.3) 

Here the trace is taken over a finite subspace of the collective Hilbert space. This 
assumption allows us to utilize the invariance of the trace under cyclic permutation, a 
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property that plays an essential role in the manipulation of the second term of ( f4.3| ). 
In this term, A is a Lagrange multiplier operator associated with a set of constraints 
that are the non- vanishing canonical commutation relations. These, at least, must 
be imposed because, to start with, the variational quantities are arbitrary matrix 
elements within the space included in the trace, and these variations must be subject 
to the kinematical constraints. The fact that we can limit the constraints to the 
commutators for canonical pairs is, at the moment, an observation. 

In order to come out with the Heisenberg equations of motion upon variation, it 
is necessary, as well, to choose variations of the coordinates that depend only on the 
coordinates themselves, so as to validate the formula 



5V(Q) = {dV{Q)/dQ k )5Q k . (4.4) 
The equations of motion that follow from these assumptions are 

-ill?, A] = (4.5) 

i[fl,A] = f?. (4.6) 

To make these agree with Heisenberg's equations, the choice A = H c r is indicated. 

We apply this formalism to the case with non-vanishing Berry potentials but con- 
stant mass matrix, described by the Hamiltonian operator, 

H ca = + V, (4.7) 

where it is understood that V may include the correction terms discussed in the previ- 
ous section. The operator equations of motion that follow from ( J4.5| ), (03), and ( f4.7| ), 
after some standard manipulations, are 

Q* = DjB ji , (4.8) 
A = -Vt + ^iQi,^}, (4.9) 



where 



and 



M = \{(dAi/dQ 3 ),Q j }i (4-10) 
= (dAi/dQ j ) — (dAj/dQ*), (4.11) 

A = Q J Bji, (4.12) 



with Bji the matrix inverse to B^ 1 , is the solution of Eq. ( [4.8| ) for the gauge- covariant 
momentum. 
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From the combination of ( [4.9| ) and ( }4.12| ), we can derive the familiar Lagrangian 
equations of motion 

(>% + V i - \{T 3h Q j } = 0. (4.13) 

Of some interest is that the associated variational principle from which these equations 
can be derived, namely 

= 5Tr(-L eS ) 

= STri-iQ^Bji-Q'Ai+V), (4.14) 



can be obtained directly from the starting trace variational principle, ( |4.3j ), simply by 
eliminating the momenta in favor of the velocities by means of the relations given above. 
Instead of independent variations of coordinates and momenta, only the coordinates 
are to be varied. In this variation, the velocity Q l is to be replaced by — i[Q\ H e $} 
and if e ff is not to be varied. One obtains the correct equations of motion because 
algebraic manipulation of the resulting trace expressions, containing commutators with 
H ef[ , that is aimed at isolating the coordinate variations, produces the same results as 
integration by parts of the time derivatives in the corresponding classical variational 
principle, proper attention being paid to the order of non-commuting factors. For the 
case just concluded, the outcome is thus completely satisfactory. 

When we reinstate the Q dependence of the mass tensor, an additional problem 
is encountered. If we insist on the cyclic invariance of the trace, which has played an 
essential role in the previous applications of the trace variational principle, then we 
have 

TriP^iP^B^}} = ■n\{r,r l .ir>}. (4.15) 

However, direct transformation, using the commutation relations, yields 

{/',. {I'r B ij }} = 2{ /', /', /,"' } + B%, (4.16) 

leading to an apparent contradiction upon formation of the trace. Of course, in deriving 
ft4.16|) , we have used the commutation relations appropriate to an infinite Hilbert space. 



The point is the same as that TrQP ^ TrPQ for a truly canonical pair. In order to 
derive the equations of motion from a trace variational principle, we deal with this 
problem by adopting ( |4.15[ ), but including the last term in ( [4.16| ) as part of the potential 



energy. This implies the assumption that the trace is taken, to start with, only over 
a finite-dimensional vector space and that therefore the invariance of the trace under 
cyclic permutation is a correct operation. As a consequence, the equations of motion 
are given in matrix form. The operator form of the equations of motion is recognized 
as a limit of these equations. In this limit, the application of ( |4.16|) after variation, can 
be seen to yield the equations of motion in the desired symmetrical (Weyl) structure. 
Thus the final variational principle takes the form 

STriH^-iHesiQ^Pi]}, (4.17) 
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where 

H mod = H cS + l -B%. (4.18) 

In the discussion above, we have noted a problem only in connection with the 
part of the kinetic energy quadratic in the momenta. It may be verified that a similar 
difficulty does not occur for the terms that are bilinear in the momenta and the "vector 
potential" . 

Unfortunately, in the transition from Hamilton's equations to Lagrange's equa- 
tions for this case additional complications are encountered; because of the non- 
commutativity of the mass tensor with the momentum, the Lagrange equations do 
not follow from the expression obtained by eliminating the momenta in favor of the ve- 
locities in the starting variational principle. The equations of motion that follow from 
this conjectured Lagrangian variational expression differ from the correct equations of 
motion by additional "quantum" potentials. This line of inquiry does not appear to 
be illuminating and therefore will not be pursued. 

5 Illustrative models 

It is simplest to illustrate the main points by choosing models in which the collective 
coordinates have already been identified, so that we need not enter into the intricacies 
of the theory of large amplitude collective motion per se. For instance, in the first 
model studied below, there are no terms in the Hamiltonian linear in the fast variables. 
Thus the model satisfies the decoupling conditions exactly. In both of the models to 
be studied, we break time-reversal invariance, and in that sense they are somewhat 
artificial. 

5.1 The first model: Berry phase in excited states 

We study the Hamiltonian 



H — H COTC + H sp + -f^i n t 

= H COTC + H N c (5.1) 

Ycore = ^ + V(Q), (5.2) 

H sp = u(Q)(a\a 1 + aia 2 + l), (5.3) 

= -G(Q-a\a 2 + Q+alai)- (5-4) 



Here the coordinates Q and the canonical momenta P are both two-dimensional vec- 
tors. We use the notation Q = (Qi,^), Q± = Qi ± iQ-2, and Q = (Q\ + Ql) 1/2 - 
(In this section, we are not maintaining the distinction between upper and lower in- 
dices.) Furthermore, the ai,a\,i = 1,2 are boson destruction and creation operators, 
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G is a coupling strength, and the frequency, u)(Q), of the uncoupled boson modes has 
been given a so far unspecified dependence on Q that will be chosen for analytic and 
numerical convenience. Note that the Hamiltonian, ( p.lp conserves the boson number, 

N = a\a\ + c4a2 = constant. (5.5) 

Since the N = problem is completely trivial, the first interesting case for the 
present model is N — 1. Here, the state vectors may be written exactly as a superpo- 
sition 

|n) = J dQ{(Q|nl)aJ|0] + (Q|n2)4|0]}, (5.6) 

where |0] is the vacuum state. The use of square brackets for the vacuum state of 
the fast variables is consistent with the notation introduced in Sec. [| The resulting 
eigenvalue equation in the space of the collective variables, Q, is determined by the 
two by two effective Hamiltonian, H e ff, with matrix elements 

(H e ff)n = H core + 2u(Q) = (H e ff) 22 , 

(H eff ) 12 = -GQ_ = (H eff y 21 . (5.7) 

Below we shall describe our solutions of the associated Schrodinger equation. 

These exact solutions are to be compared with the adiabatic approximation, with 
and without the Berry potential terms. For this approximation, we require the normal 
modes of H^c, which we calculate in a standard way from the equations of motion, 

[ai,H NC \ = ujai-GQ_a 2 , 

[a 2 ,H NC ] = ua 2 -GQ + a 1 , (5.8) 
by forming the matrix elements, 

A = [OH*]. (5.9) 

With Q representing the energy of the state |^], we obtain the equations 

Qipx = uj^ x - GQ_ip 2 , 

= -GQ+fa+wfo, (5.10) 



that yield the eigenvalues 



ft«(Q) = u(Q)-GQ, 

Q {2 \Q) = u(Q) + GQ, (5.11) 



that are degenerate when Q — 0. The associated normalized solutions of (|5.10|) are 
represented most conveniently by introducing the normal-mode creation operators 

61 = ^4 (5.12) 
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In detail we have possible choice, 

b\ = -j=a[ + -^=expz0(Q)4, (5.13) 

h\ = -L exp ^(Q)a t 1 --^=4, (5.14) 

where 

tan0(Q) = Q2/Q1. (5.15) 

We now apply these elementary results to the adiabatic approximation. In this case 
we represent a suitable subset of the eigenfunctions ( |5.6|) in the form 

|n) = J dQ(Q|n)6}|0:Q], (5.16) 

where the notation : Q refers to the vacuum for the normal modes. (For the current 
model, it coincides with the uncoupled vacuum.) The considerations of Sec. |3| now 
apply to this class of state vectors and, in particular, we apply Eq. ( |3.23|) . As a special 
case of this equation, we have 

i i i 

where 

Ai = {[O-.Qlb&bllO-.Q], (5.18) 

(4)21 = i[0:Q\b2dib{\0:Q], (5.19) 



where <9j means the partial derivative with respect to Qi. With the help of Eqs. ( |5.13| ) 
and ( |5.14| ), the quantities of interest are found to take the values 



A = ^(g 2 ,-Qi), (5.20) 
El(^) 2 i| 2 = (1/4Q 2 )- (5.21) 

i 

The singular character of these results was to be expected. 

The collective or adiabatic Hamiltonian, He, that thus emerges from the assump- 
tion that the state vectors of interest can be written in the form (|5.16|) , has the structure 



H c = (1/2M)[(P - A) 2 + (A 21 ) 2 ] + V(Q)) + (3/2)fi«(Q) + (1/2)^(Q). (5.22) 

We turn to the problem of solving the associated eigenvalue problem. We wish 
to compare the results of exact (numerical) calculations with the eigenvalues of the 
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collective Hamiltonian. It is useful to choose u(Q) such that the latter is exactly 
solvable. One such choice is 

uj(Q) = l -GQ. (5.23) 
In this case the additional contribution to the potential is zero, 

+ V)(Q) = 0. (5.24) 

We further take M = 1 and V(Q) = \Q 2 . 

We use the definition of the two-dimensional angular momentum, 

T - 9 

= -i(QiPi-Q2Pi), (5.25) 

to simplify the P • A term. We can further simplify the equation H c ijj(Q) = Etp(Q) by 
substituting 

HQ) = Q~ 1/2 x(Q)e im * (5.26) 

and find 

HI? + w im2 + m + j 1 + \ Q2 } xiQ) = ExiQ) (527) 

Even though the centrifugal term has changed, this is still very similar to the radial 
equation for the two-dimensional harmonic oscillator. It can be solved by the substi- 
tution 

x{q) = q a+1,2 z- q2l2 L^\ q 2 ). (5.28) 
As is well known the right hand side of this equation satisfies the condition 

d 2 1 - 4a 2 

~d^ x + (g + ~w~ )x = (4n + 2a + 2)x ' (5 ' 29) 

We thus find that 
and 



a{m) = ^1/4+ (m + 1/2) 2 (5.30) 



E nm = (2n + a(m) + 1). (5.31) 
Without the Berry's phase terms we would have found 

E nm = (2n+ \m\ + 1). (5.32) 

When G = the solution ( p.32[ ) is exact. Thus Eq. ( |5.31| ), which is independent of G, 
cannot be valid for all G. It should be valid in the adiabatic limit, which means that 
the two frequencies Q must be very different. This occurs when G is very large. 
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The exact solution of the problem can be calculated using a spherical harmonic 
oscillator basis for the Q coordinates. This is coupled to a state containing either one 
a\ or one a 2 boson, 

\n,m,n 1 ,n 2 ) = \n, m)(a t 1 ) (ni) (4) {n2) |0]. (5.33) 

The interaction Hamiltonian only couples states with n\ = l,m — mi to states with 
n 2 = l,m — mi + 1. If we thus choose the value mi we have only a relatively small 
matrix to diagonalize. Since we wish to obtain accurate results for large G, we allow 
the number of Q harmonic oscillator quanta to be fairly large. We have used up to 200 
harmonic oscillator states in each block (which leads to a 400 by 400 matrix eigenvalue 
problem). 

In table [I] we give a selected set of results for m\ — and a number of values of G. 
Similar results for mi=l are listed in table [2| We clearly see the convergence to the 
the collective model results including the Berry phase for large G. 

Table 1: The lowest 10 eigenvalues of the coupled problem for mi = as a function 
of G. The column labeled Berry lists the eigen energies of the collective Hamiltonian 
02j 



G 


0.01 


0.1 


1 


10 


100 


1000 


Berry 


Ei 


1.00875 


1.07849 


1.36229 


1.58463 


1.66532 


1.69302 


1.7071 


E 2 


2.01328 


2.13134 


2.92007 


3.48844 


3.63505 


3.68298 


3.7071 


E 3 


3.01540 


3.14492 


4.09426 


5.39473 


5.60870 


5.67438 


5.7071 


E 4 


4.01827 


4.18125 


5.30157 


7.29413 


7.58443 


7.66658 


7.7071 


E 5 


5.01997 


5.19061 


6.57000 


9.17666 


9.56145 


9.65932 


9.7071 


E 6 


6.02221 


6.22068 


7.66964 


11.0262 


11.5393 


11.6525 


11.7071 


E 7 


7.02367 


7.22768 


8.90918 


12.8133 


13.5177 


13.6459 


13.7071 


E 8 


8.02558 


8.25425 


10.0230 


14.4981 


15.4965 


15.6396 


15.7071 


E 9 


9.02687 


9.25972 


11.1883 


16.0926 


17.4754 


17.6335 


17.7071 


Eio 


10.0286 


10.2840 


12.3410 


17.7119 


19.4544 


19.6276 


19.7071 



5.2 The second model: ground state 

We next study a model that is also exactly decoupled, differing from the one just 
investigated only in the form of the interaction. In this model Eqs. fl5.1|) -( |5l3[ ) stand 
as given, but Eq. ( |5.4j ) is modified to 

Hint = -G Q(a\a 2 + a\ai) - Gi(Q + aia 2 + Q-a\a\). (5.34) 
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Table 2: The lowest 10 eigenvalues of the coupled problem for mi = 1 as a function 
of G. The column labeled Berry lists the eigen energies of the collective Hamiltonian 
ET22. 



G 


0.01 


0.1 


1 


10 


100 


1000 


Berry 


Ei 


2.01308 


2.11333 


2.41566 


2.54921 


2.57601 


2.58042 


2.58114 


E 2 


3.01670 


3.17369 


4.20360 


4.51926 


4.57057 


4.57959 


4.58114 


E 3 


4.01807 


4.16340 


5.44213 


6.48766 


6.56479 


6.57866 


6.58114 




5.02086 


5.21512 


6.43129 


8.45380 


8.55880 


8.57766 


8.58114 


Er, 


6.02202 


6.20307 


7.86356 


10.4165 


10.5527 


10.5766 


10.58114 


E 6 


7.02436 


7.25000 


8.79958 


12.3739 


12.5464 


12.5755 


12.58114 


E 7 


8.02538 


8.23691 


10.1070 


14.3218 


14.5400 


14.5744 


14.58114 


Eg 


9.02745 


9.28063 


11.1975 


16.2512 


16.5336 


16.5733 


16.58114 


Eg 


10.0284 


10.2669 


12.3141 


18.1346 


18.5270 


18.5721 


18.58114 


Eio 


11.0302 


11.3082 


13.5364 


19.8627 


20.5204 


20.5709 


20.58114 



The second term of this interaction spoils the conservation of boson number used to 
simplify the solution of the previous problem. For the present problem, it is already 
interesting to study the spectrum when the fast variables are in their ground levels, 
since there are now non-trivial ground-state correlations and associated non-trivial 
Berry potentials, that we shall calculate. 

We first turn to the study of the BO approximation. As before we need the normal 
modes of the fast variables, as determined from the the equations of motion 



[d, H NC ] = ooai — G Qa 2 — GiQ^a 



[a 2 ,H NC ] = ua 2 — GqQcl! - GiQ^a[, (5.35) 
and their hermitian conjugates. In terms of the definitions 

b\ = tfjfa^-xfa,, (5.36) 



i>f = [0:QK-|*»], (5.37) 
X f = [0:Q|o}|*W], (5.38) 
= b\\0:Q], (5.39) 



we obtain the eigenvalue equations 

flip! = uipi - GoQip 2 - 



fi^2 = ^2 - G Qipi - GiQ^Xh 
-fiXi = ^xi - G Qx2 - GiQ + ipi, 

-Q X 2 = ux 3 -G Qxi-G 1 Q + if> 1 . (5.40) 
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The physical eigenvalues are the positive roots of 



n 2 (Q) = (00 t G QY - GjQ\ (5.41) 

which are again degenerate at Q = 0. As for the first model studied, we label these so- 
lutions fiW, i = 1,2. The corresponding amplitudes are determined from the equations 
of motion and the normalization conditions 

hM 2 + |^| 2 -|Xi| 2 -|x 2 | 2 = l. (5.42) 

The simplest forms for the amplitudes are achieved by repeated use of the eigenvalue 
equation ( |5.41| ). We thus find 

4 j) = (-iy' + Vi j) =real, (5.43) 
xF = exp(i<J))x?, (5.44) 
xf = (-iy +1 x? =re a \, (5.45) 

V1 ^[{u + ^y + i-^G^ + n^ + ^Gl-G^Q 2 }^ V ' ; 
X? = (-1)'" +1 / C, ' ) (WP ) , (5.47) 
fCi) = 9lQ (5 48) 

Here is the polar angle for Q, <fi = aictaji(Q 2 /Qi). 

These results are thus available for application to the BO approximation studied 
by means of the assumption 

|n) = J dQ|Q)(Q|n)|0:Q]. (5.49) 

It follows that the Berry vector potential is given by the formula 

A i = i[U'.Q\d i \U'.Q}. (5.50) 

To carry out this calculation, we need the form of the correlated vacuum state, as given 
by the equation [ ]H) | 

\0:Q]=Mex P [^Z ll a]a}}\0}, (5.51) 

where M is a real normalization factor, whose value we shall not need, Oj|0] = 0, and 
Zij is a generally complex array whose values will be considered below. Using the 
constancy of the norm of the correlated ground state as a function of the values of Q, 
we can rewrite (|5.50|) as 



A = U{[0:Q\di\0:Q] - [8,0 : Q\0 : Q}} , (5.52) 
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which shows that the (real) normalization factor M does not contribute. This leaves 
the result 

A = -i[(a\a\)diZ k i - {a k ai)diZ^, (5.53) 

where the averages, indicated by angular brackets, are with respect to the correlated 
vacuum. 

We next indicate how (|5.53|) can be evaluated in terms of the solutions we have 
found above for the equations of motion. The quantities Z k i are solutions of the 
equations 

that follow from the definition 6j|0:Q] = 0. The expectation values can also be evalu- 
ated with the help of the formulas 



^?Z* k = x[\ (5.54) 



o i = V?' ) 6 i + X? ) *6j, (5-55) 

and their hermitian conjugates, that are the inverse to Eq. (|5.36|) . An elementary 
calculation now yields 

(a[a\) = (a kai y = (5.56) 
Thus for the vacuum Berry potential, we obtain the formula 

A = ^[x k m) ri m> d t ZM - xV^t^iZU- (5-57) 

Any other matrix element, (Ai) u0 of the Berry potential can be calculated by similar 
techniques. In fact, the only non-vanishing elements of this type are (Aj)n i0 , (^1)20,0, 
and (Aj)o2,o, where, for example, v = 20 means a state with two correlated bosons of 
type 1 and none of type 2. We quote formulas for these matrix elements that follow 
from the same elementary techniques that furnished Eq. (|5.57|) . We thus have 

(A)n,o = i[U:Q\di\0:Q] 
= -i[di 11:Q|0:Q] 
= -z[0:Q|Mi(6 2 )|0:Q]} 

= -i[Y.4 >d ^f*-xf*d^f\ (5.58) 
j 

Similarly, remembering the normalization of the states, we find 

(AKo = ^[E^xf* - xf^d^f*]. (5.59) 

Finally, the matrix element (A)o2,o is obtained from 

(|5.59| ) simply by replacing all superscripts 1 by 2. In order to implement fully the 
BO approximation, it is necessary to evaluate these matrix elements in terms of the 
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explicit solutions of the model displayed in Eqs. ( |5.43| )-( [5^46| ). It is also appropriate to 
remark here that in consequence of these relations and Eq. (|5.54j) we can write 

Zj k = exp(—i(j))Zj k , Zj k = real (5.60) 

The calculation of the Berry potentials as defined above is now straightforward. In 
terms of purely real quantities, we find, for example, 

A = (d i( j>) [(x?¥x } + X?¥i ] )Zn + (x { M 1} ~ xfVP^ia]. 
Using the real form of Eq. ( |5.54| ), this expression can be simplified to 

G\ 



(5.61) 



+ 



1 

(u + nW)(u + fi« - 2G Q) + (Gg - G\)Q 2 
1 



+ + fl(2) + 2G Q) + {Gl - G?)Q 2 



:-Q2,q 



It- 



(5.62) 



In contrast to the previous model the vector potential is not singular at the origin, but 
take the finite value 

G\ 



A(0) 



;(-Q 2 ,Qi). 



Turning to the off-diagonal matrix elements, we find by similar calculations 



(-4)n,o 
(4)20,0 



0. 



(5.63) 



(5.64) 



v / 2exp(-z0)[(^0)4 1) X^ 1) + K^diX? ~ xPwP] 
-v / 2exp(-^) (rf 1} ) 2 / (1) [d l <p + i(d l Q)d Q \nf^' 
^2™-' ,^ G.m^ + lu) - G Q] 



exp^— u, 



2[(u + + - 2G Q) + {Gl - G?)Q 2 ] 



x 



^(-Q 2 ,Qi) + i-jj(Qi,Q2)Qd Q lnfU 



(5.65) 



(A 



i 02,0 



-V2 



exp(— i(f 



G^GoQ + in^ 



2[(u + + - 2G Q) + {Gl - G\)Q 2 } 



x 



^(-Q*,Qi) +i^(QuQ2)Qd Q \nf^ 



(5.66) 



These off-diagonal potentials are also regular as Q goes to zero, for reasonably well 
behaved uj{Q). 

We now study the numerical solution of the adiabatic Hamiltonian, 
H C = h(P - A) 2 + |A 20 , | 2 + |A 02 ,o| 2 ] + \Q 2 + \P 1 \Q) + W\Q)], (5.67) 
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2 4 6 8 10 

Q 

Figure 1: Parameters in the collective Hamiltonian for G = 1, ujq = 2. The lower 
panel shows the frequencies tt^, the middle panel shows the size of the diagonal Berry 
potential, and the upper panel shows the square of the off-diagonal Berry-potentials. 

where we have once again made the choice M = 1,V(Q) = \Q 2 - For simplicity 
we take Go — 0, and write G\ = G. For this special choice = = Q and 
(^4;)o2,o — — (A)20,o- A simple form for u(Q), chosen such that Q is positive definite, 
is 

cu(Q) = u + GQ. (5.68) 

In Fig. [I] we show some of the relevant quantities in the collective Hamiltonian for 
the choice G — 1, Uq — 2. We diagonalize the collective Hamiltonian by first going 
to spherical coordinates, and use the fact that m is conserved to write down a radial 
equation. We then map the Q values from the interval [0, oo) to the interval [0,1]. 
Finally we make a finite difference approximation to the radial equation, and solve the 
approximate equation by matrix diagonalization. 

Solving the complete problem, without making the adiabatic approximation is 



24 



somewhat involved. Since the Hamiltonian is invariant under the interchange a\ <-> a 2 , 
we introduce new operators that are invariant under this parity transformation, 

ai = ^=(a\±4). (5.69) 

The interaction Hamiltonian takes the simple form 

H int = -G„(4a + - ala_) - ^(Q+[a 2 - a 2 ] + Q-[^ ~ at 2 ]). (5.70) 
We now make the expansion 

, 2K-2n / + \ 2k 

a! 



(4V , ■ 

* = E^(g)-S= == lo] (5.7i) 



x£ ^(2tf-2«)!(2«)! 

(The even powers in this equation constitute the only combination that includes the 
vacuum for the fast degrees of freedom.) 

An approximate solution is obtained by limiting the sum over K, while summing 
over all allowed k. At the same time we expand iPk,k{Q) in a finite number of spherical 
harmonic oscillator eigenfunctions (Q\nm). We denote the m value used in the expan- 
sion of i/jk,k(Q) by m K,K- The interaction does not couple states of different k for fixed 
K. We thus find that m^ K = = mx-i — 1. Thus the value m is a constant of 
the motion, and can be used to specify different solutions. This quantity corresponds 
directly to the value of m in the collective Hamiltonian. We have performed matrix 
diagonalizations for ujq = 2, 10 and G — 1. We have used harmonic oscillator states 
up to principal quantum number 70, and K values up to 20. (This corresponds to a 
61,401 by 61,401 matrix.) The resulting matrix, which is very sparse, was diagonalized 
using a Lanczos algorithm. The eigenvalues were checked for convergence by compar- 
ing a calculation with smaller cut-offs on n and K against one with larger cutoffs. In 
tables |3] and § we compare a few selected ground state energies of the complete collec- 
tive Hamiltonian with the exact solution. The splitting between states with opposite 
values of m is completely due to the vector potential. For the case u> = 2, where 
the adiabatic approximation we see that the size of the splitting is very close to the 
exact value. The difference is probably due to non-adiabatic effects that can not be 
completely neglected for uj = 2. For the case ujq = 10 the correspondence is much 
closer, but the size of the Berry vector potential is much smaller as well. 
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Table 3: A comparison between the exact numerical ground state energies and those 
for the collective Hamiltonian. u = 2,G = 1, and the values of m are listed in the 
table. 





m = —5 


m = 5 


m — — 1 


m = 1 


exact 


9.54715 


9.73101 


4.97309 


5.03405 


collective 


9.56722 


9.74656 


4.98669 


5.04375 



Table 4: A comparison between the exact numerical ground state energies and those 
for the collective Hamiltonian. ujq — 10, G — 1, and the values of m are listed in the 
table. 





m 


= -5 


m 


= 5 


exact 


18 


.0745 


18 


.0918 


collective 


18 


.0750 


18 


.0922 



A Scalar product and transformation of momen- 
tum operators 



In this appendix we provide the proof of Eq. ( [2.7|) , which specifies how the momen- 
tum operator transforms under a general point transformation. It will be essential to 
recognize that this result is tied to a choice of scalar product. We suppose that the 
Hamiltonian ( |2.1| ) is to be used in conjunction with the metric 

(* a |*»> = Jde---deK(0M0- (A.i) 

Thus 

ir a - -i{d/d£ a ). (A.2) 
Now carry out the point transformation ( |2.3| ) with Jacobian J, 

J=\dC/dx l3 \. (A.3) 

If we introduce a new wave function 

= jh a , (A.4) 

the metric is preserved in the sense 

(*„,* b ) = (^aiA) = [ dx 1 ---dx N r a {x)^h{x). (A.5) 
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This is the metric that is associated with Eq. ( |2.7| ), as we proceed to show. 
With d^ 1 ■ ■ ■ d£ N = [<%], we study 



[d£]*>a*6 = J [dx\r a J^aJ~^b- (A.6) 



Thus we must show that 



We have first 



j\ir a r\ = \{f (x1 {-id/dx i )}. (A.7) 
ill Qx l Id l 

jw* = j- 5 (a.8) 

For ease of notation, consider below a three-dimensional space, since this will exhibit 
the general features. By utilizing the equation, equivalent to the definition of J as a 
determinant, 

with 6ijk the alternating symbol, we display the algebraic manipulation leading to the 
desired result, 

i i 1 __I 1 dC d __I 

J2n a j 2 = -J 2 -e ijk e a p 7 — J 2 

l OX J OX K ox 1 

119 l dtfd? 1 

1 1 9 x d^d? 1 
" 2T 2 W^ l ^ 7 ^a? j 2j 



1 _i, d£P dn d 

(A.10) 

By an application of (|A.9|), this becomes 



I I 9 , I (9x* . life' 9 I 1 Ji ! 9 , . . 

27 2 &^ 2 &p + 2i 2 d^'dx i 2 ~ Ytd^'dx 1 ^ ^ ' ^ 



the result sought. 



B Potential energy for the fast variables 

In terms of the original coordinates, let us consider the change in the potential energy 
between two neighboring points, £ and £ + 5£. To second order in 5£ we have the usual 
terms of a Taylor expansion, 

AV = v{Z + 6Z)-V{Z) 

= vjc + lv^csf, (B.l) 
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that now appears most unsatisfactory, since AV is a scalar, but the second term of 
(B.l| ) contains the ordinary rather than the covariant second derivative. This defect is 
removed by replacing the quantity 5£ a by substitution from the relation 



(B.2) 



that contains the Christoffel symbol, 



r 



Pi - 



\ B a& (-8,5/3,7 + B$7,P ~ -8/37,-5) • 



(B.3) 



We thus obtain the form 



AV = V a dC + lV a0 dCdf, 



(B.4) 



where the second term contains the covariant derivative 




transformation to any alternative set of coordinates such as the is standard. 

This allows us to calculate the quantity of Eq. (2.10) from given dynamical 
quantities, provided we can define a complete set of coordinate axes at each point of 
the decoupled manifold. This is done in two stages. In the first stage, the basis vectors 
at each point of the tangent space to the collective submanifold, for example the set 
fa, are determined by the algorithm that discovers the collective submanifold. A set 
of basis vectors orthogonal to the tangent space is then determined (non-uniquely) 
by the requirement that these be orthogonal to the f l a and to each other with respect 
to the metric B a @. The basis vectors orthogonal to the tangent space are precisely the 
elements needed to compute the Hamiltonian of the fast variables. 
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